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ABSTRACT 

The effects of pinning between fiuxoids and vortices in the core of a neutron 
star, on the dynamics of the core neutron superfluid are considered. The pinning 
impedes, but does not absolutely block, any radial as well as azimuthal motion 
of the neutron vortices with respect to the lattice of fiuxoids. The time scale for 
the coupling of rotation of the core superfluid to the rest of the star is calculated, 
allowing for the effect of the finite frictional force on the neutron vortices due to 
their pinning with the fiuxoids. This turns out to be the dominant mechanism 
for the coupling of the core of a neutron star to its crust, as compared to the 
role of electron scattering, for most cases of interest. Furthermore, different 
behaviors for the post-glitch response of the core superfluid are distinguished 
that might be tested against the relevant observational data. Also, a conceptually 
important case (and controversial too, in the earlier studies on the role of the 
crustal superfluid) is realized where a superfluid may remain decoupled in spite 
of an spinning up of its vortices. 

Subject headings: stars: neutron - hydrodynamics - pulsars 



1. INTRODUCTION 

The neutron superfluid in the core of a neutron star coexists with the super-conducting 
protons and the "normal" degenerate electrons. The interior charged plasma including the 
lattice of the proton vortices (fiuxoids) is expected to be strongly coupled to the lattice 
of nuclei and the electrons i n the crust due to the strong magnetic field present, with a 



coupling time scale < 10 s (lAlpar et al.l Il984bl ). All charged components of the star are 



hence considered as one co-rotating "component" , referred to as the "crust" ; in contrast to 
the neutron superfluid in the core, being the second "component" . An isolated neutron star is 
subject to an electromagnetic spin-down torque that acts on its magnetic field, hence on the 
"crust" . The core superfluid would be likewise driven to follow the long-term spinning down 
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of the crust through the mutual coupling mechanism that would operate between the crust 
and the superfluid vortices. The force, on the vortices, responsible for such a spinning down 
(or any assumed short-term relaxation) of the core superfluid has been, generally, considered 
to be only that of the scattering of the electrons off the vortex cores. The dominant scattering 
effect is that of the induced magnetization of the vortices caused by a "drag" between the 
proton and neutron condensates. The associated vorte x velocity relaxatio n time scale is 
r e ~ 10P S - 20P S , where P s is the spin period of the star (lAlpar et al.lll984al ). 



However, in the quantum liquid interior of a neutron star a neutron vortex is expected 
to "pin" to a fluxoid should the two structures overlap. The mechanism of the pinning is 
associated with either the proton density perturbations or the magnetized nature of both 
the fluxoids as well as the neutron vortices. The strength Ep of the corresponding energy 
barrier has been estimated as Ep ~ 0.1 or 1.0 MeV, for the two mechanisms, respectively. 
Likewise, the effective length scale dp of the pinning interaction is expected to be of the order 
of the coherence length £ p or the London penetration depth A p of the proton superconductor 
for the above two mechanisms, respectively. The effective pinning force (per intersection) 
fp = is therefore roughly the same for the two interaction mechanisms, since the value 
of Ep due to the magnetic interaction is larger than that of the densit y perturbation by the 
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The consequences of such a pinning with respect to the radial motion of the vor- 
tices have been discussed in the literature, to some extent. The outward moving vor- 
tices in a slowing down neutron star have been argued to sweep t he fluxoids along with 



them and hence expel the magnetic flux out of the stellar core (ISrinivasan et al.l 11990 



Jones! Il991bl ; ISedrakian fc Sedrakianl Il992l ) . The magnetic evolution of pulsars, particu 



larly those in binary systems, has been calculated, based on this effect ( Ding et al 



1993 



Jahan-Miri & Bhattacharvalll994j ; [jahan-MirilE996l : IUrpin et all997l ; [jahan-Mirill200oh . Also, 
due to the steady-state slowing down of the star, the superfluid would be rotating faster than 
the vortices, resulting in a radial Magnus "force" that drives them against the pinning bar- 
riers. Glitch inducing mechanisms, caused by the core sup erfluid, have been suggest ed on 
this ground, invoking the possibility of sudden disturbances (IMuslimov &; Tsyganlll985l ). and 
based on the predicted long-term evolution of the st eady-state value of the associated ro- 
tational lag between the superfluid and the vortices (I Jahan-Miril |2002| ) . Also, the expected 
post-glitch behavior of the obs ervable rotation f requency of the crust of a neutron star has 
been discussed in this context (IChau et al.lll992l ). 



The pinning would however also impede any assumed relative azimuthal motion between 
the two lattices of the vortex lines, and impart azimuthal forces. That is, it could impart 
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a torque, and act as a means for a transfer of angular momentum to the core superfluid. 
The effect is, in principle, same as the role of the standard "static" frictional force, realized 
i n the laboratory experimen t s on a rotating sup erfluid, in presence of the pinning centers 
( ITsakadze Tsakadzd Il980l ; lAdams et al.l Il985f ) . Such a role for th e azimuthal comp onent 
of the pinning force in the c ore of a neutron s_tar_was referred to by iJahan-Miril ( 119981 ) , and 
has been recently studied by lSidery fc Alparl (120091 ). These authors assume however absolute 
pinning, and assign an infinite strength to the azimuthal component of the pinning force, 
all owing o nly fo r a sliding motion of the vortices along the floxoids, as was first suggested 
by ISaulsl (119891 ) . It should be clarified, at the outset, that we are concerned with only 
an application of the standard results a nd basic for mulation of the superfluid rotational 
dynamics, at the hydrodynamical level (jSoninl 119871 ). avoiding complications due to the 
microscopic details. We adopt the generally accepted standard structure for the interior of a 
neutron star. We attempt to include the role of the azimuthal force on the vortices, due to 
their pinning with the fluxoids, in the rotational dynamics of the core superfluid, in addition 
to the viscous drag force, due to the scattering of the electrons. 

In the following we discuss the rotational dynamics of neutron stars taking into account 
also the above effects (both radial, as well as azimuthal) due to the pinning of the vortices 
in the core of the star. The steady-state relative rotation of the different components of a 
neutron star is briefly described, distinguishing between the cases of pinned and unpinned 
vortices in the core superfluid. It may be noted, that a clear picture of the steady-state 
relative rotations is indeed needed in order to follow the discussion further. The vortex- 
velocity relaxation time scale is then calculated, allowing also for the finite effect of the 
frictional force of the pinning with the fluxoids. The dynamical timescale for the coupling 
of the pinned core superfluid is thence derived. In addition, qualitatively different behaviors 
for the response of the core superfluid to a jump in the rotation frequency of the "crust" are 
distinguished, depending on the magnitude of the jump, in the case of pinned superfluid. In 
the final section, the relevance of our results to the observational data on pulsars glitches 
are discussed briefly, indicating the possibility of distinguishing among the different cases 
considered. 



The CORE SUPERFLUID 



Superfluid vortices move with the local superfluid velocity except when there is an 
external force acting on a vortex. For a given ext ernal force F ex , per unit length of the 
vortex, the equation of motion is given (ISoninlll987l ) as 



p s K x [v s - v L ) 
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where v s and vi are velocities of the local superfluid and the vortex line, respectively, p s is 
the superfluid mass density, k is the vorticity of the vortex line directed along the rotation 
axis, with a magnitude k = for the neutron superfluid, where m n is the mass of a 
neutron. The kinematic side of the equation is generally referred to as the Magnus "force". 
The force, arising from the gradient of the superfluid kinetic energy, is exerted between the 
superfluid and the vortices, and explains how a "mutual" friction between the superfluid 
and its environment is realized. As is usual in this field, we will neglect the collective effects 
of the vortex lattice, finite tension of a vortex, as well as any relevant non-dissipative force 
on a vortex b y the lattic e, which has been studied recently for the vortices in the crust of 



neutron stars (jLinkll2009l ) 



In the case of a rotating superfluid, the number density n v , per unit area, of the vortices 
obeys nn v = 2Q S , and the rate Cl s of change of the rotation frequency Q s of the superfluid is 
associated with a radial velocity v r of the vortices, (l s = —2 — v r , where r is the distance 
from the rotation axis, and v r > is in the outward direction. Thus, a spinning down (up) 
torque on a rotating superfluid is accompanied by a radial outward (inward) motion of the 
vortices. 



2.1. Steady-State Relative Rotations 

During the steady state spinning down of a neutron star, the vortices are expected to 
be co-rotating with the crust, including the proton condensate and the lattice of fluxoids, in 
the core of the star. Strict co-rotation of the vortices with the crust is not however possible 
for an assumed vortex relaxation process due only to the electron scattering. In contrast, the 
pinning force could be imparted even while the vortices are co-rotating with the fluxoids. The 
pinning potential barrier could act as a "stretched spring" , that does pull the mass attached 
to it, even while both are moving together. In any case, the difference in the rotation rate 
between the vortices and the fluxoids (the crust) is tiny and will be neglected, for the steady 
state considered. On the other hand, the steady-state rotation rate of the superfluid should 
be different than that of the crust and/or the vortices, due to the assumed pinning of the 
vortices. The superfluid rotational lag is defined as u> = Q s — Ql, where is the rotation 
frequency of the neutron vortex lines. 
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2.1.1. The Steady-state Lag 

The superfluid critical lag u cr is determined from a balance of the radial Magnus force 
Fyir — PsKtu with the radial component of the exter nal force, F r , on t he vortices, per unit 



length. Equating F r with the pinning force results in (jDing et al.lll993[ ) 



W cr = ~~~~T ( 2 ) 



p s K,rdf 



B 1 ' 2 



~ 1.6 x lO-^rads- 1 ) =^-, (3) 

where r& is the value of r in units of 10 6 cm, and Byi is the strength of the magnetic field B c 
in the stellar core in units of 10 12 G. For the typical magnetic fields of young pulsars, one 
finds typical values of 

i) w cr > 10~ 4 rad s~\ 

for i? 6 ~ 1, at the outer regions of the stellar core, and 

H) w cr > 1CT 3 rad s _1 , 

for Rq ~ 0.1, at the inner parts, embracing ~ 1% of the moment of inertia of the core. 



c • 



The steady-state spinning down will be thus established while Q c = < Q s , Cl s = Cl 
and the steady-state vale of the lag = uj ct > 0, with the above estimates for its value 
(which also indicates a state of differential rotation for the core superfluid at the different 
radial distances). It may be noted that random unpinning of the vortices could in principle 
serve to establish even smaller values of the lag < u CT . An estimate of the effect, base d 



on the vortex creep model, shows however the difference to be negligible (IChau et al.lll992l ). 
Thus, we will be neglecting the role of random unpinning of the vortices, and hence assuming 
Woo = u cr , in the following. 

The expected relative rotation of the "crust", the core superfluid, and its vortices is 
sketched in Fig. 1, as a visual aid which turns out to be vital too. The sketch helps to 
keep track and be ensured that the, preliminary but, fundamental physics applicable to a 
transfer of angular momentum between the superfluid and its environment (the crust, here) 
is obeyed! That is, for a superfluid spin-down (-up) to be achieved the crust must be, or 
tend to be, rotating slower (faster) than not only the superfluid itself, but also slower (faster) 
than the vortices (see §3.2.2 and §4, below). Each plot indicates the three successive stages 
at a disturbance and the transitions between them. Namely, the conditions during a steady 
state, as discussed above (marked as "S.S."), followed by a jump in the rotation frequency 
of the crust ("JUMP"), and after the core has responded to that jump, returning to a new 
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steady state ("post-G"). The cases of pinned and free vortices behave differently and are 
shown separately, as there is also a further distinction for the pinned cases depending on the 
relative size of the jump, to be discussed below. 



2.2. Vortex Velocity Relaxation 

An assumed departure from the steady-state co-rotation of the vortices with the fluxoids 
brings about their continually crossing one another, with an energy cost of the pinning energy 
Ep, per intersection, associated to a frictional force fp. The total number Ax of crossings 
per relative rotation cycle of the two lattices (embedded in a spherical boundary surface) is 
iVx = N{N V sin% (see the Appendix), where JVf and N v are the total number of the fluxoids 
and the vortices, and x is the angle between the two families of lines, namely the angle of 
inclination between the rotation and the magnetic axes of the star. Excluding the cases 
of parallel lattices corresponding to \ ~ or 180 deg which are not common among the 
observed pulsars, the correction due to the sin% factor may be neglected, assuming near 
perpendicular geometry between the two lattices, for simplicity. 

The effective frictional force F P [ n of the pinning, per unit length of a vortex, while cutting 
through the lattice of uniformly distributed fluxoids (with a spacing df) may be estimated 
by taking a time average over a cycle ( "pinning cycle" ) of passing through a pinning region, 
of a length li = dp, and an adjacent inter-pinning zone, of a length l 2 = d { — dp. This is in 
contrast to the drag force due to the electron scattering F e i ec , which acts continuously, and 
will be treated separately. Thus, 

F pin = -^-(i 1 / 1 +* 2 / 2 ) (4) 

J cyc 

where T cyc = ti + t 2 is the total time for a vortex to pass through one pinning cycle, and the 
subscripts refer to the two regions, respectively. The pinning force is operative only within 
each pinning zone, namely 

fi = N c f P (5) 

= - = ^, (6) 
d{ dp df 

h = 0, (7) 

where we have used N c = j- for the number N c of the potential crossings per unit length of 
a vortex, per pinning cycle (see the Appendix). 



Next, the microscopic motion of any given vortex, while cutting through the fluxoids, 



-7- 



might be expected to fall in between the following two limiting types of behaviors, for a first 
approximation. 



1. A vortex might preserve the same (relative) velocity (with the ffuxoids) in both the 
above two regions, throughout a pinning cycle. Hence, the travel times, for the two 
regions, would be linearly proportional to their sizes. Namely, 



t\ l\ dj 



t 2 h df — dp' 
ti dp 



T C yc d{ 



and (8) 

(9) 



which results in 



^Pin = ^ (10) 

~ 6 x 10 12 (dyn cm" 1 ) E McY B l2 (11) 
where E MeV is the value of E P in units of MeV. 

2. Or, the vortex might be expected to adjust its velocity instantaneously as it moves 
through the pinning and the free regions, successively. The vortices might therefore 
spend almost all of the time inside the pinning regions and fly across the free spacings 
between the fluxoids rapidly, in almost zero time. That is 

t 2 = 0, and (12) 
ti = T cyc (13) 



which results in 



~ 1.4 x 10 16 (dyn cnT 1 ) £ McV d^i (15) 

12 



where d\ 2 is the value of dp in units of 10 cm. 



The above two cases might each represent a better approximation for either a transient 
relaxation, following a sudden change of the rotation frequency of the fluxoids (as for a 
glitch) or during an equilibrium steady state spin-down phase, respectively. In either case, 
the force is independent of the relative velocity of the two families of lines; a dynamical 
feature different than that of, say, the force due to the electron scattering. 
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The corresponding vortex relaxation time Tp, needed for a given relative velocity be- 
tween the fluxoids and vortices to be dissipated, may be defined, from the EOM of a unit 
volume of the charged component gas co- rotating with the fluxoids (the "crust"), through 

B v~t-v c 

n v F pin = pc—; (16) 

-tp 

where v c is the local velocity of the electron gas and the fluxoids, ie. the "crust", and p c 
is its effective density, corresponding to the contribution of the real crust together with the 
free protons (and the electrons) in the moment of inertia of the star. Note that since F P j n is 
a velocity-independent force, T P is the total time for the decay of the velocity, and not an 
exponential time constant. For an assumed jump AQ C in the rotation frequency Q c of the 
crust, out of its steady state co-rotation with the vortices, one derives, from Eq. 16 together 
with Eq. 11 or 15, 

T P = 830 (s) A ^ c - Pl " , or (17) 



0.4 



Att c p 13 d 12 



for the two approximations considered above, respectively, and where pis is the value of p c 
in units of 10 13 gem -3 . Notice that the dependence on the distance r from the rotation axis 
(v = rfl) has been averaged out, and the radial velocity of the vortices, during a relaxation, 
would have negligible effect. 

Substituting typical values of ~ 10~ 6 , and E MeV ~ p 13 ~ B 12 ~ 1, one finds, 
from Eq. 17, Tp ~ 10 -4 s, which is muc h shorter than the corresponding time scale due 



to the electron scattering r e ~ 1 - 2 s flAlpar et al.l Il984bl ). The estimated value of Tp 



for the other case, from Eq. 18, is even smaller, by more than three orders of magnitudes. 
Nevertheless, because T P does depend on the initial value of the induced relative velocity 
it may, in principle, become larger than r e , but that might happen only for the very large 
assumed disturbances, much larger than that observed even in the giant glitches of the Vela 
pulsar. For the steady state spin down, on the other hand, Tp would be even much smaller 
than the above estimate. 

The drag force, T elec , of the electron scattering off neutron vortices is, likewise, associated 
with a vortex relaxation time r e , through an equation similar to Eq. 16, as is the relation 
between the overall relaxation time scale r v with the total force T to t = T P i n + F e i ec , namely 

v~i-v c 

n v F tot = p c , (19) 



which means 



111 

~ = 7fT + - (20) 
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Therefore, the frictional pinning force turns out to be the dominant coupling mechanism for 
the core superfluid, and t v ~ T P may be adopted, in general. 



3. Coupling Timescale of the Superfluid 



The dynamical timescale for the rotational coupling of the bulk superfluid to its envi- 
ronment (the crust) is the time needed for the simultaneous re-adjustment of the vortices 
in both radial and azimuthal directions in response to an external torque on the superfluid, 
exerted primarily on the vortices. The external fo rces on the yortice s could, in ge neral, 
be of a viscous drag or a "static" frictional nature (jAdams et al.lll985t lJoneslll991bl ). For 
the superfluid in the core of a neutron star, both types are present and correspond to the 
electron scattering and the fluxoid "scattering", discussed above, respectively. The latter 
type, associated with the "pinning" forces, should not be however confused with the role of 
the pinning forces on the pinned vortices co-rotating with the pinning centers. In order for 
the pinning forces to act as frictional forces and impart a net torque on the superfluid the 
vortices should maintain a radial velocity, hen ce unpinning continuously due to the effect 
of the Magnus force, which requires \u\ > u cr (jAdams et ajj [1985| ; I Jahan-Miril l2005a| ) . For 



| a; | < u> cr , random unpinning events might play a role in the long-term coupling of the su- 
perfluid, but the effect is neglected here, for simplicity and more so because of its negligible 
effects on the quick post-glitch responses. A more general treatment should include the 
super fluid coupling rate driven also by the random unpinning of the vortices (IJahan-Miri 
20061 ). as has be en already studied, in the context of the vortex creep model, and found to 
be unimportant (jSidery fc Alparl 120091 ) . 



3.1. The case of free vortices (no pinning barriers) 



The dynamical timescale Tf ree of the superfluid rotational relaxati on, in the absence of 



any pinning of th e vortices (Fig. la), has been previously determined (lAfpar fc Saulslll988 
Jahan-Miril ll998T ) as 



r /2p2 
_ fs -'c-'s 

Tfree — t e ' ia IT J 

1 lbn z Ll 



(21) 



This is the exponential timescale that appears in the solutions obtained for the relaxations 
of the radial r v (t) and the azimuthal <fi v (t) components of the vortex positions, in polar 
coordinates on the equatorial plane, as a function of time t. The time behavior of the vortex 
position (following an assumed sudden rise in the rotation frequency of the "crust" , say at a 
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glitch) is governed by the vortex equation of motion (Eq. 1), using F ex = F e iec, in the absence 
of any pinning. The solutions, which were not given in the correct form previously, are: 



r v (t) 



To 



sf) 



o + tt F t + 



+ 1- 



n 



sf) 



-t/Tfr, 



1/2 



r, - 



Ikti v 



In 



r 



(22) 
(23) 



where 0-subscripts indicate initial values at t — of the corresponding quantities, / = I c + I s , 
and Q F = j(I c Q c0 — I s ^so) is the final equilibrium frequency of both the crust and the 
superfluid. 



3.2. The case of Pinned Superfluid 



For a pinned superfluid, already in a steady-state with fi c = i\ < ^s and co = = 
u CT > 0, as discussed in § 2.1, a sudden jump in Q c is then assumed to occur, the cause 
of which is irrelevant for the purpose of the present discussion. The crossing through the 
fluxoids by the vortices is inevitable until the state of co-rotation between them is reached 
(considering rigid vortex lines), hence the vortices would be subject to forward azimuthal 
pinning forces. The response of the pinned superfluid to such disturbances is argued, below, 
to be qualitatively different, depending on the magnitude of the increase AQ C in the rotation 
frequency of the crust. The effect, should not be however attributed to the particular form(s) 
of the pinning force, calculated above. Indeed, the same types of behaviors would be expected 
even in the presence of the drag force due to the electron scattering alone, as long as the 
neutron vortices are assumed to be subject to the pinning to the "crust" (co-rotating with 
the fluxoids and the charged components). Also, note that the decomposition of the effect 
into a "jump" in Q c and a "subsequent response" of the superfluid is not meant to be in real 
time, rather it is for the sake of the analysis, as is common. 



3.2.1. Large Jump 

Following an assumed "large" jump, such that AQ c0 > u CT (Fig. lc), the fluxoids (and 
the electron gas) would be rotating faster than both the vortices and the superfluid itself. 
Hence, the existing forward azimuthal forces on the vortices would be consistent with a 
corresponding inward motion of the vortices required for a spin-up torque to be imparted to 
the superfluid. The steady state lag would be thus washed out during (the jump and) the 
following quick relaxation of the superfluid (and the vortices), bringing the superfluid, first, 
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to a state of co- rotation with the fluxoids and the vortices. Hence, following a large jump 
the superfluid will be, first, spun up by the crust, on a dynamical time scale r pin , until the 
equilibrium state Q s = = Q c is reached (see Fig. lc). 

The coupling timescale T p - m may be determined from a calculation similar to the above 
case of the free vortices, except that r v (Eq. 20) would now be the relevant vortex relaxation 
time. That is (compare with Eq. 21) 

I I 2 P 2 

Tp in = - r v + -^- r v -\ (24) 
pm I v 16tt 2 IJ v V ; 

It may be noted that, the pinning force, fp, at each intersecting point, is directed solely along 
the direction perpendic ular to the intersecting fluxoid; t he component of the force along 



the fluxoid being zero (jSaulsl Il989l ; ISidery &: Alparl I2009T ) . In solving the vortex equation 



of motion (Eq. 1), with the substitution F ex = F pin , a decomposition of the force to the 
radial and azimuthal components would be initially faced with a difficulty, that may be 
however bypassed. The projections of the fluxoids, on the equatorial plane, form an array of 
parallel lines. In the presence of a difference between the rotation rates of the vortices and 
the fluxoids, the vortices at any given radial distance r cut through those lines at varying 
angles. Hence, the azimuthal component of the pinning force, at each intersection, would 
vary between its maximum vale fp down to zero, for the vortices located at the different 
azimuth angles. The same would be true for the radial component of the force, as well. Thus, 
at any given time and for the vortices at any given radial distance, the average values of the 
azimuthal and the radial components of the pinning force, per intersection, may be used; 
averaged over the vortices at varying azimuthal angles over a circle. The average values, 
would be the same for both components and, simply amount to - fp. This will correspond 
to a change in the relaxation time Tp by a factor |, that we have neglected. 



3.2.2. Small Jump 

In contrast, following an assumed "small" jump, in the rotation rate of the crust, such 
that AQ c0 < u cr (Fig. lb), there would follow no transfer of angular momentum between the 
crust and the superfluid, initiated by the jump. The superfluid could be neither spun down 
nor up, by the crust! Following the jump, the vortices would be indeed spun up to come 
into a co-rotation with the fluxoids, due to the existing forces. Nevertheless, the superfluid 
which is already rotating faster than the crust, and its own vortices, (Fig. lb) could not be 
possibly spun up, and gain angular momentum from the slower component, during such a 
spin-up of the vortices. Also, the superfluid could not be spun down, during such a spinning 
up of the crust, that is accompanied by a spinning up of the vortices. Because, a spin-down 
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of the superfluid requires a radially outward motion of the vortices. Nonetheless, an outward 
motion of the vortices may be realized only if they do rotate, or tend to be rotating, faster 
than the t he superfluid envir onment, so that a slowing down torque could be imparted on 
them (see Jahan-Miri J2005al lbh for a more extended explanation; indeed false conclusions, 



about the behavior of the superfluid in the crust of a neutron star, have been communicated 
because of the neglect of this primitive but fundamental fact, as discussed therein). 

Hence, following a small jump, the co-rotation of the vortices with the crust is achieved 
with no transfer of angular momentum to the superfluid, ie. at constant Q s . The time scale 
involved would be only that of the azimuthal vortex velocity relaxation, Tp. The speeding 
up of the vortices alone does not, however, require any transfer of angular momentum, since 
these are but massless (super)fluid configurations (neglecting the inertia of the vortices, as 
is usually assumed). 



4. Discussion and Observational Implications 



Glitches are observed in radio pulsars as sudden changes AQ C in the rotation frequency 
f2 c of the crust with observed values of the jump in the range 10 -9 < ^ 10~ 6 . In younger 
pulsars, the jump in Q c is also accompanied by an incre ase Afl r in t he observed spin-dow n 



rate fl c of the crust, with typical values of < 2.5% rtLynd 119871 : iKrawczyk et al.l 120031 1 . 



The effect, as such, has been explained successfully in terms of a decoupling of the superfluid 
i n the crust of a neu tron star, which has a fractional moment of inertia of similar magnitudes 
(lAlpar et al.lll984al ). However, there has been observed cases with 



> 10% 



(25) 



and recovery timescales up to ~ 44 d (jLynd 119871 ; iFlanaganl Il995l ) . These observed cases 
would necessarily imply that part of the star with the same fractional moment of inertia( ie. 
> 10%, and indeed up to 60% at some observations) has been decoupled from the crust, at 
the time of the observation. This, by itself, is a definite proof of a rotational decoupling of (a 
part of) the core of the star, at those glitches, since the crust does not constitute that much 
of the moment of inertia of the star. Given also the observed recovery time scales at these 
events, the effect, in turn, indicates that a pinne d superfluid compon ent of the stellar core 
must be involved, as has been argued previously (j Jahan-Miril l2005bl ) . The core superfluid, 



in the absence of pinning, does similarly cause observed large values of Ap, but that would 
last only over a time of the order of Tf ree (Eq. 21). In contrast, the pinned superfluid remains 
decoupled until the superfluid rotational lag, which is washed out quickly at the jump (see 
Figs lb & lc), recovers its steady state value u CT . During this recovery period tpc, the 
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observable crust is being spun down, by the unchanged external torque N ext , at an increased 
post-glitch rate ^pg, that could be as large as, 

fipG = f^ss, (26) 

where Cl^s = ^f 1 is the steady state value of the spin down rate of the star. Since, the 
vortices are pinned and co-rotating with the fluxoids and the crust, thus the superfluid is 
decoupled and the lag is built up due only to the resulting decrease in f2 c . Hence, the 
post-glitch recovery timescale tpc, over which the increased spin-down rate of the crust may 
persist, is expected to be 

tpc = -^ >10 5 (s), (27) 

"PG 

where the estimated value is for an assumed low value of u cr = 10~ 4 rads _1 , a large value 
of j- = 20, and a typical value of (l c = 10~ 10 rads~ 2 . Obviously, larger recovery times are 
expected for the lower values of fipc, that would be expected if only a fraction of the core 
superfluid takes part in the decoupling-coupling mechanism. Therefore, the largest observed 
values for the post-glitch spin-down rates of the pulsars, as well as the associated recovery 
timescales recorded, find a natural explanation in terms of the expected decoupling of the 
superfluid component in the core, being subject to the pinning with the fluxoids. 

Further, the two cases of the large versus the small jumps, considered above, would have 
different implications for the observable post-glitch recovery. In the case of a large jump, the 
expected initial spin-up of the core superfluid, discussed above, would show up as a further 
increase, by an additional factor j-, in the observable spin-down rate of the crust, compared 
to that given in Eq. 26. In contrast, no such an increase is expected for the small jumps 
(compare Fig. lb with lc). Such a difference in behavior, between the two cases, would 
however persist only over time periods ~ r pin , which makes it indeed hard to be detected, 
observationally. 

On the other hand, the two cases might be judged and distinguished based also on the 
earliest detected values of Af2 c , as compared to the value of u cr . That is, the two cases 
of small versus large jumps correspond to (initial) values of AQ C < u cr and AQ C > u cr , 
respectively (as indicated in Figs lb and lc). Notice that, for a large jump the condition 
AQ C > uj cr would be still true, even after the initial relaxation of the superfluid, over the 
time Tpi n , ie. at the post-glitch t = on Fig. lc. Also, the observational uncertainty in the 
exact epoch of the glitch, being of the order of few minutes, would not blur between the two 
cases, since the change in Q c over such periods would be negligible. Therefore, given the 
generally observed initial values for 4pS cited above, as compared to the expected values of 
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w cr ^ 10 4 rad s 1 for the pinned core superfluid, one may conclude that the case of a small 
jump is the relevant one for most, if not all, of the observed glitches. 

Moreover, the two cases would have different implications also for the possible cause of 
the glitches. A large jump has to be necessarily induced by some "external" agent, other than 
the core superfluid itself, simply because a donor of the angular momentum could not spin up 
its counterpart to a frequency more than its own. In contrast, and by the same token, a glitch 
induced by the core superfluid would necessarily rank as a small one. The core superfluid 
is thus, statistically, more favored as the potential cause of the small jumps. The above 
mentioned indication of the data for the prevalence of the small jump s, would thus further 



support the glitch indu cing mechanisms driven by the core superfluid ( IMuslimov fc Tsygan 



19851 : 1 Jahan-Miril 120021 ) . 



Finally, the pinned core superfluid does not respond to a jump in Q c , at a glitch, smaller 
than the critical value of the superfluid rotational lag, u CT , as indicated above. The angular 
momentum gained by the crust, at such glitches, would not be shared, any further, with the 
core superfluid. The case is highlighted here again, not because of any prominent observable 
effect different than the other case of the large jumps. Rather, for its prominent conceptual 
value in demonstrating the vital requirement for the presence of a corresponding torque for 
a spin-down (-up) of a superfluid to be achieved. The motion of the superfluid vortices too 
has to conform with this basic physics, be it a smooth motion in the absence of any pinning, 
or the so-called creeping of the pinned vortices. 

This work was supported by a grant from the Research Committee of Shiraz University. 



A. The average rate of crossing with fiuxoids, per unit length of a vortex 

We are assuming the fiuxoids, in the core of a neutron star, form a uniform array of 
parallel lines, along the magnetic axis (B) of the star, having a total number of iVf, within a 
spherical boundary surface of radius R. Likewise, the neutron vortices constitute a uniform 
lattice of parallel lines, along the rotation axis (fl) of the star, having a total number of N v , 
within the same spherical boundary. The two axes are in general inclined at an angle Xi °f 
which we consider first the case with x — f > f° r further clarity, and will generalize at the 
end. Taking the z-axis along Q, the fiuxoids are thus parallel to the x-y plane. 

1. A circular slab, perpendicular to the z-axis, its rim at the polar angle 9, a distance 
z = RcosO from the origin, with a radius R z = RsinO, and a thickness dz = sin 9 d9, 
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includes a number N\> of fluxoids, being parallel to it, where 

= - N f sin 2 6 d9 (A2) 

2. The fraction of the fluxoids in the slab, lying at (cylindrical) distances r to r + dr from 
the z-axis is 

N a (r) = ^-N h (A3) 

3. Each of the lines at r, sweeps an annular area s(r), in the x-y plane, per rotation cycle, 
that is 

s(r) = 2(tiR 2 - nr 2 ) (A4) 

4. The area dS h swept by all the fluxoids in the slab would be 

dS h = 2 f RZ N s (r)s(r) (A5) 
Jo 

= ^nR z 2 N h (A6) 

5. The cumulative total area St t swept by all the fluxoids within the sphere, parallel to 
the x-y plane, thus becomes 

Stot = 2 / 2 dS * ( A7 ) 
= 7rR 2 N{ (A8) 

This may be compared with a corresponding value ^nR 2 N{ which would result for an 
assumed cylindrical geometry, or with a value 2ttR 2 N{ if all the fluxoids were of the 
same length 2R, lying radially. 

6. The number Ax of crossings between the fluxoids and the vortices (the latter all passing 
through the x-y plane, with a number density n v = ^j), per cycle of relative rotation, 
per unit area of the x-y plane, becomes 

A^x = n v S tot = AT v 7V f (A9) 

which may be, in turn, compared with a corresponding value 2N V N{, if each fluxoid 
were to cross each and any of the vortices, per relative rotation cycle. 
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7. The total sum L v of the lengths of all the vortices, within the sphere, is 

L v = n v (^7r^ 3 ) = ^RN V (AlO) 

Hence, the average number Aq of crossings, per unit length of a vortex, per rotation 
cycle, is 

8. For an square lattice of the fluxoids, with a spacing df, between the nearest fluxoids, 
Thus, 



"• = -3- (A12) 



* = iT (A13) 



9. The point of intersection of each vortex with the x-y plane, at a distance r from the 
z-axis, moves a distance d{r) = 2irr, per rotation cycle. The corresponding average 
distance d v that a vortex moves, becomes 



1 



R 



d v = - d(r) dr = nR (A14) 
R 

10. Therefore, the rate N c of crossings, per unit length of a vortex, per pinning cycle, is 
derived as 

N c = ^N 1 = 3 -± (A15) 
d v 4 df 

11. For the general case of an arbitrary inclination angle x, the derivation would be similar, 
except that one starts by considering a slab, parallel to the fluxoids and at the same 
angle x with respect to the rotation axis, so that it preserves its angle during the 
rotation. Then the projected area swept by those fluxoids, on the x-y plane would 
have an additional factor sinx, as compared to that in Eq. A8. As a consequence, 
the total number of the crossings would be corrected by the same factor, namely 
Ax = A" v N { sinx, and also for the final result, A c = | ^ sinx- 
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Fig. 1. — Schematic representation of the relative values of the angular velocities of the 
superfluid Q s (dotted line), the vortices I^l [dashed line), and the "crust" Q c (full line), 
excluding the overall spinning down of the system. On each plot, three successive phases are 
shown at a relaxation, namely an earlier steady state (marked as "S.S."), an assumed sudden 
increase in f2 c ("JUMP") and the subsequent relaxed state ("POST-G"). a) is for the case 
of free vortices, in the absence of pinning barriers, whence the steady state spinning down 
corresponds to (an almost) co-rotation of the vortices, with the superfluid and the crust. 
The assumed increase in the angular velocity of the crust is shared with the core superfluid 
over its dynamical timescale td = Tf ree . b) and c) are both for a pinned superfluid. Due to 
the assumed pinning with flux lines, in the initial spinning-down steady-state, a rotational 
lag between the superfluid and the vortices is necessary. For the "small" jumps, as in b), 
only the vortices are spun-up over an "azimuthal" relaxation timescale Tp^ = Tp. For the 
"large" jumps, as in c), the pinned core superfluid couples to the crust over its corresponding 
dynamical timescale td = r pin , as discussed in the text. 




time — > 



